o 

(N 

< 

m 



00 



c^ 



Symmetry, Integrability and Geometry: Methods and Applications SIGIVIA 6 (2010), 070, 17 pages 

C-Integrability Test for Discrete Equations 
via Multiple Scale Expansions 

Christian SCIMITERNA and Decio LEVI 

Dipartimento di Ingegneria Elettronica, Universita degli Studi Roma Tre and Sezione INFN, 
Roma Tre, Via della Vasca Navale 84, OOI46 Roma, Raly 



E-mail: scimitema@fis. uniromaS. it , levi@roma3. infn. it 



Received May 29, 2010, in final form August 20, 2010; Published online August 31, 2010 



doi:10.3842/SIGMA.2010.070 



Abstract. In this paper we are extending the well known integrability theorems obtained by 
multiple scale techniques to the case of linearizable difference equations. As an example we 
apply the theory to the case of a differential-difference dispersive equation of the Burgers 
hierarchy which via a discrete Hopf-Cole transformation reduces to a linear differential 
difference equation. In this case the equation satisfies the Ai, A2 and A3 linearizability 
conditions. We then consider its discretization. To get a dispersive equation we substitute 
the time derivative by its symmetric discretization. When we apply to this nonlinear partial 
difference equation the multiple scale expansion we find out that the lowest order non- 
secularity condition is given by a non-integrable nonlinear Schrodinger equation. Thus 
showing that this discretized Burgers equation is neither linearizable not integrable. 
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• ' Calogero in 1991 [3] introduced the notion of S and C integrable equations to denote those 

^— V . nonlinear partial differential equations which are solvable through an inverse Spectral transform 

(^ I or linearizable through a Change of variables. Using the multiple scale reductive technique he 

was able to show that the Nonlinear Schrodinger Equation (NLSE) 

- idtf + d,,f = p2\f\^f, f = f{x,t), (1.1) 

^ ' appears as a universal equation governing the evolution of slowly varying packets of quasi- 

monochromatic waves in weakly nonlinear media featuring dispersion. Calogero and Eckhaus 
then showed that a necessary condition for the S'-integrability of a dispersive nonlinear partial 
differential equation is that its multiple scale expansion around a slowly varying packet of 
quasi-monochromatic wave should provide at the lowest order in the perturbation parameter an 
integrable NLSE. Then, by a paradox, they showed that a C-integrable equation must reduce 
to a linear equation or another C-integrable equation as the Eckhaus equation [H [22] . 

In the case of discrete equations it has been shown [26l[Il[l6l[l8l[I71[l9l[ini[IIl[l2] that 
a similar situation is also true. One presents the equivalent of the Calogero-Eckhaus theorem 
stating that a necessary condition for a nonlinear dispersive partial difference equation to be 
5-integrable is that the lowest order multiple scale expansion on C'°°) functions give rise to 
integrable NLSE. The nonlinear dispersive partial difference equation will be C-integrable if 
its multiple scale expansion on C^°°' functions will give rise at the lowest order to a linear or 
a C-integrable differential equation. 

By going to a higher order in the expansion, the multiple scale techniques give more stringent 
conditions which have been used to find new 5-integrable Partial Differential Equations (PDE's) 
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and to prove the integrability of new nonlinear equations [8l[9l[T5]. Probably the most important 
example of such nonlinear PDE is the Degasperis-Procesi equation [7]. Up to our knowledge 
higher order expansions for nonlinear linearizable equations have not been considered in details. 

The purpose of this paper is to show that the integrability theorem, stated in [11], can 
be extended to the case of linearizable difference equations, providing a way to discriminate 
between S-integrable, C-integrable and non-integrable lattice equations. The continuous (and 
thus discrete) higher order C-integrability conditions are, up to our knowledge, presented here for 
the first time. We apply here the resulting linearizability conditions to a differential-difference 
dispersive nonlinear equation of the discrete Burgers hierarchy [21] and its difference-difference 
analogue. 

In Section [2] we present the differential-difference linearizable nonlinear dispersive Burgers 
and its partial difference analogue and discuss the tools necessary to carry out the multiple scale 
C-integrability test. In Section [3] we apply them to the two equations previously introduced, 
leaving to the Appendix all details of the calculation, and present in Section |4] some conclusive 
remarks. 

2 Multiple scale perturbation reduction of Burgers equations 

The Burgers equation, the simplest nonlinear equation for the study of gas dynamics with heat 
conduction and viscous effect, was introduced by Burgers in 1948 [2l[27j. Explicit solutions of 
the Cauchy problem on the infinite line for the Burgers equation may be obtained by the Hopf- 
Cole transform, introduced independently by Hopf and Cole in 1950 [UIH]. This transformation 
linearize the equation and the solution of the linearized equation provide solutions of the Burgers 
equation. 

Bruschi, Levi and Ragnisco [21] extended the Hopf~Cole transformation to construct hier- 
archies of linearizable nonlinear matrix PDE's, nonlinear differential-difference equations and 
difference-difference equations. 

The simplest differential-difference nonlinear dispersive equation of the Burgers hierarchy is 
given by 

dtUn (0 = ^ {[1 + h^n (t)] [Un+l {t) " U„ (t)] - ^^^^i^-^j , (2.1) 

where the function n„(t) and the lattice parameter h are all supposed to be real. Equation ()2.ip 
has a nonlinear dispersion relation u = _£Hi^|l^_ 

When /i ^- and n — ?> oo in such a way that x = nh is finite, the Burgers equation ()2.ip 
reduces to the one dimensional wave equation dtu — dxU = 0{h^). 

Through the discrete Cole~Hopf transformation 

hcpn (t) 

equation ()2.ip linearizes to the discrete linear wave equation 

Mn [t) = ^ . (2.3) 



The transformation ()2.2p can be inverted and gives 

j=n-l 






Yl (l + huj), n>no + l, (2.4a) 
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^"0 



j=no-l 

n (i+H) 

j=n 



n < uq — 1, 



(2.4b) 



where (pno = (pnoit) is the function 0„ calculated at a given initial point n = uq. When n„(i) 
satisfies equation ()2.ip . the function </>n(i) will satisfy the discrete wave equation ()2.3p if (/)„(, 
satisfies the ordinary differential equation 



1 
2h 



1 + hun 



1 



(1 + hUn-l] 



0, 



n=no 



whose solution is given by 

(pnoit) =(?i„o(*o)exp|— / 



1 + hUn 



1 



(1 + hUn~l] 



dt' 



n=no 



n— >— oo 



and to is an initial time. If lim n„(t) = n„oo is finite equations ()2.4p reduce to 

1 + hu,.. 



7=n— 1 

cj)n = ait){l + hU-ocT n 

7=— oo 

where a (t) is a t-dependent function, 



1 + hu. 



«(t) = aoexp|l(l + ^^i— - 
A C-integrable discretization of equation ()2.ip is given by the partial difference equation 



^n,m- 



-1 -Itr^ 



O" 



1 

2h 



(1 + hUn,m){Un+l,m — ""n.m.+l) 



^71— l,m ^n,m+l 



1 + hUn- 



l.m 



(2.5) 



where a is the constant lattice parameter in the time variable. Equation (j2.5p is dissipative 
as it has a complex dispersion relation uj = Mn [l + i^ sin (k/i)] . As we are not able to con- 
struct a dispersive counterpart of equation (|2.ip we consider a straightforward discretization of 
equation p.ip 



2^ 



1 
2/1 



(1 + hUn,m) {Un+l,m 



[Un~l,m Un^m) 
1 + hUn-l,m 



(2.6) 



whose nonlinear dispersion relation is 
sin (k/i) a 



sin (ojcr) 



/i 



In the remaining part of this section we will present the tools necessary to carry out the 
multiple scale expansion of equations ()2.ip , ()2.6p and construct the conditions which they must 
satisfy to be C-integrable equation of j order, with j = 1, 2, 3, i.e. such that asymptotically they 
reduce to a linear equation up to terms respectively of the third, fourth and fifth order in the 
perturbation paramether. These conditions up to our knowledge have been presented for the 
first time by Dr. Scimiterna in his PhD Thesis [25] and are published here for the first time. 

2.1 Expansion of real dispersive partial difference equations 

For completeness we briefly illustrate here all the ingredients of the reductive perturbative 
technique necessary to treat difference equations, as presented in [lOl [TTl [25] . 
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2.1.1 Prom shifts to derivatives 

Let us consider a function «„ : Z — t- R depending on an index n € Z and let us suppose that: 

• The dependence of Un on n is reahzed through the slow variable ni = en € M, e € I 
< e <C 1, that is to say u„ = u(ni). 

• The function u (ni) € C^°°> {T>), where D € M is a region containing the point rii. 

Under these hypotheses one can write the action of the shift operator T„ such that TnUn 
Un+i = u{ni + e) as the fohowing (formal) series 

T„n(ni) = n(ni) + eu^mini) + -r-U2ni(f^i) H h -r^u^knAni) H 



= ^yu^knAni), (2.7) 

where u^kniini) = d u{ni)/dn^ = (i„^u(ni), being d^ the total derivative operator. The last 
expression suggests the following formal expansion for the differential operator T„: 



hOO 1, 



fc=0 

valid only when the series in equation (j2.7p is converging. So we must require that the radius 
of convergence of the series starting at rii is wide enough to include as an inner point at least 
the point ni + e. 

Let us introduce more complicated dependencies of Un on n. For example one can assume 
a simultaneous dependence on the fast variable n and on the slow variable ni, i.e. Un = u{n, ni). 
The action of the total shift operator T„ will now be given by r„ti„ = n„+i = u{n + l,ni + e) 

so that we can write T„ = 7n ini , where the partial shift operators Tn and Tnl are defined 
respectively by 

oo ^ 

Tji^^u{n, ni) = u{n + 1, ni) = ^ —d^u{n, ni) = e^"u{n, ni), 

fc=o 

and 

7;^')'u(n,ni) =n(n,ni + e) = ^ — (9^jn(n, ni) = e'^"in(n, ni). (2.8) 

A:=0 

The dependence of n.„ on n can be easily extended to the case of one fast variable n and K slow 
variables Uj = CjU, Cj &M, 1 < j < K each of them being defined by its own parameter Cj. The 

action of the total shift operator Tn will now be given in terms of the partial shifts 7n , Tnj, 

as In ^ In 11 Irij 

Let us now consider a nonlinear partial difference equation 



j=(-^(-),^(+)) 
l-""+'='™+iJfc=(-Ar(-),Ar(+)) 



0, (ArW,/cW)>0, (2.9) 



for a function n„^m : Z^ ^> M which now depends on two indexes n and ttt, E Z which we will 
term respectively as discrete space and iime indices. Equation (|2.9|) contains m and n-shifts, 
respectively in the intervals (m — JC^~\ m + K.^^') and {n—J\f^~>,n+J\f^'^>). Under some obvious 
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Table 1. The operators An , Bm and Cn/m appearing in equations (|2.10p 



U) 





J = 


i = i 


i = 2 


J =3 


i = 4 


Aii' 


1 


Nldn, 


2 ""1 


6 ""1 


nIq4 
24 ""1 


sH' 


1 


Ml 9™, 


^9^1 + ^29^2 


+MiM2 9™i9™2 + 
+M39„3 


24 CjTii 1 2 C'miC'm2 + 




1 


^i^'+S« 


+iViA/i9„i9„i 


+NiM2dn-,dm^ + 


+^9i,9„, + ^9^,9,„,+ 
+iViAfiM29„,9,„,9,„2 + ^^dl^d^^ + 



hypothesis on the C'°°' property of the function Un.m and on the radius of convergence of its 
Taylor expansion for all shifts in the indices n and m involved in the difference equation (j2.9p . 



we can write a series representation of Un-\-k,m+j around Un,m- We choose the slow variables as 
im m with 



^rik 



Nke\ l<k<Kn, 



M,e\ 



l<j<Kr, 



where the various constants N/^, Mj and e are all real numbers. In this presentation we assume 
Kn = 1 and K^ = K (eventually K = +00) so that 



+00 
rr _ 7-(i)7-(^"i) _ 7-(i) >p J Aii) 

i=0 

K +00 

rf _ 7-(i) TT 7-^'™.- -> _ -7-(i) V J »(j) 

/C +00 

-i-n-i-m — In 'm '"•! J_J_ 'mj — In 'm / ^ "^ ^n,mi 

j=l j=0 



(2.10a) 
(2.10b) 
(2.10c) 



where the operators An , Bfi , and Cn,m are given in Table [TJ Inserting the explicit expres- 
sions ()2.10p of the shift operators into equation (12. 9p . this turns into a PDE of infinite order. So 
we will assume for the function Un,m = u{n, m, ni, {rrij}--^ , e) a double expansion in harmonics 
and in the perturbative parameter e 



+00 7 



Un 



7=1 6l=-7 



(2.11) 



with n^ {ni,mj,j > 1) 



M 



ni , nij , j > 1) , where by a bar we denote the complex conjugate, 



in order to ensure the reality of Un,m- The index 7 is chosen > 1 so that the nonlinear terms 
of equation ()2.9p enter as a perturbation in the multiple scale expansion. For simplicity we will 

set A'^i = Mj = 1, j > 1 . Moreover we will assume that the functions u\ satisfy the asymptotic 

conditions lim uiy = 0, V7 and 6 to provide a ineaningful expansion. 

2.1.2 From derivatives to shifts 

The multiple scale approach discussed above reduces a given partial difference equation into 
a partial differential equation for the amplitudes u\ contained in the definition (j2.1ip . 
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We can rewrite the so obtained partial differential equation as a partial difference equa- 
tion inverting the expansion of the partial shift operator in term of partial derivatives (|2.8|) . 
From (|2.8p we have 

dn. = ^InTi^) = ^ln(l + 6AW) =5]L^[A(-[f , (2.12) 

fc=i 

where A„j = {Tnl — l)/e is the first forward difference operator with respect to the slow- 
variable ni. This is just one of the possible inversion formulae for the operator Tni ■ For example 
an expression similar to equation (|2.12p can be written for the first backward difference operator 

Ani = (l— [7^1 ] )/e- For the first symmetric difference operator A„-|^^ = (7^i — [7^i ] )/2e 
we get 



a„,=sinh-eAW = f ^^[aW 



k 

k=l 

where Pkifi) is the k-th. Legendre polynomial evaluated at x = 0. 

Only when we impose that the function Un is a slow-varying function of order I in the 
variable ni, i.e. that Af^^Un = 0, we can see that the dm operator, which is given by formal 
series containing in general infinite powers of the A„^ , reduces to polynomial of order at most i. 
In [TTj, choosing i = 2 for the indexes ni and mi and £ = 1 for ?7t,2, it was shown that the 
integrable lattice potential KdV equation [20] reduces to a completely discrete and local nonlinear 
Schrodinger equation which has been proved to be not integrable by singularity confinement and 
algebraic entropy |131l23j. Consequently, if one passes from derivatives to shifts, one ends up in 
general with a nonlocal partial difference equation in the slow variables n^ and m^. 

2.2 The orders beyond the Schrodinger equation 
and the C-integrabihty conditions 

The multiple scale expansion of an equation of the Burgers hierarchy on functions of infinite 
order will thus give rise to PDE's. So a multiple scale integrability test will require that a dis- 
persive equation like equation ()2.ip is C-integrable if its multiple scale expansion will go into 
the hierarchy of the Schrodinger equation 

idtip + dli) = 0. 

To be able to verify the C-integrability we need to consider in principle all the orders beyond 
the Schrodinger equation. This in general will not be possible but already a few orders beyond 
the Schrodinger equation might be sufficient to verify if the equation is linearizable or not. In 
the case of S-integrable nonlinear PDE's the first attempt to go beyond the NLSE order has 
been presented by Degasperis, Manakov and Santini in [8]. These authors, starting from an 
S-integrable model, through a combination of an asymptotic functional analysis and spectral 
methods, succeeded in removing all the secular terms from the reduced equations order by order. 
Their results could be summarized as follows: 

1. The number of slow-time variables required for the amplitudes u- appearing in (|2.1ip 
coincides with the number of nonvanishing coefficients of the Taylor expansion of the 
dispersion relation, ujj {k) = 4 ^j}j ' . 

2. The amplitude u\^ evolves at the slow-times rUg, s > 2 according to the s-th equation of 
the NLS hierarchy. 
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3. The amplitudes of the higher perturbations of the first harmonic u^ , j >2 evolve, taking 
into account some asymptotic boundary conditions, at the slow-times m^, s > 2 according 
to certain linear, nonhomogeneous equations. 

Then they concluded that the cancellation at each stage of the perturbation process of all 
the secular terms is a sufficient condition to uniquely fix the evolution equations followed by 
every Uj\ j >1 for each slow-time rus ■ Point 2 implies that a hierarchy of integrable equations 
provide for a function u always compatible evolutions, i.e. the equations in its hierarchy are 
generalized symmetries of each other. In this way this procedure provides necessary and sufficient 
conditions to get secularity-free reduced equations [8]. 

We apply the present procedure to the case of C-integrable partial difference equations. 
Following Degasperis and Procesi [9] we state the following theorem: 



Theorem 1. // equation (12. 9p is C -integrable then, after a multiple scale expansion, the func- 
tions Uj , j > 1 satisfy the equations 

d„,J^^ - i-iy-^Bsd^^y,'^ = M,nW = 0, (2.13a) 

Ujf = fS), (2.13b) 

Vj, s > 2, where Bgd'^^u^ is the s-th flow in the linear Schrodinger hierarchy and Eg are real 
constants. All the other u^ , k>2 are expressed as differential monomials of Ur , r < j — 1. 

In equation ()2.13bp fs{j) is a nonhomogeneous nonlinear forcing written in term of differential 
monomials oi Ur , r < j. From Theorem [T] it follows that a nonlinear partial difference equation 
is said to be C-integrable if its asymptotic multiple scale expansion is given by a uniform 
asymptotic series whose leading harmonic u^^' possesses an infinity of generalized symmetries 
evolving at different times and given by commuting linear equations. Equations (I2.13P are 
a necessary condition for C-integrability. 

It is worthwhile to stress here the non completely obvious fact that, in contrast to the first 
order wave equation, dfU = dxU, all the symmetries of the Schrodinger equation commuting 
with it are given by the equations (I2.13ap and only by them. This implies that all the equations 
appearing in the multiple scale expansion for a C-integrable equation are uniquely defined. 

It is obvious that the operators Mg defined in equation (I2.13ap commute among themselves. 
However the compatibility of equations ()2.13bp is not always guaranteed but is subject to some 
compatibility conditions among their r.h.s. terms fs{j)- Once we fix the index j > 2 in the set 
of equations (|2.13bp . this commutativity condition implies the compatibility conditions 

Msfs' (j) = M,,f, (j) , Vs, s' > 2, (2.14) 

where, as fs {j) and fs' (j) are functions of the different perturbations Ur of the fundamental 
harmonic up to degree j — 1, the time derivatives dm^ > dm^, appearing respectively in M^ and Ms' 
have to be eliminated using the evolution equations (j2.13p up to the index j — 1. These last 
commutativity conditions turn out to be a linearizability test. 

Following [8] we conjecture that the relations (|2.13p are a sufficient condition for the C- 
integrability or that the C-integrability is a necessary condition to have a multiple scale expan- 
sion where equations ()2.13p are satisfied. To characterize the functions /s(j) we introduce the 
following definitions: 

Definition 1. A differential monomial M [u, ] , j > 1 in the functions u, , its complex conju- 
gate and its ni-derivatives is a monomial of "gauge" 1 if it possesses the transformation property 

M [uf^] = e'M {uf\ , when uf = e^'uf . 
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Definition 2. A finite dimensional vector space Vy, v >2 \s the set of all differential polyno- 
mials of gauge 1 in the functions u- , j > 1, their complex conjugates and their ni-derivatives 
such that their total order in e is v, i.e. 



Older (^d^^u- ) = order (5^^ u^- ) = ^ + j = v, /i > 0. 

Definition 3. Vy{ii), fi > 1 and z^ > 2 is the subspace of Vu whose elements are differential 
polynomials of gauge 1 in the functions u, , their complex conjugates and their ni-derivatives 
such that their total order is u and 1 < i < /i. 

From Definition [3] it follows that Vu = Vu{i'—'2)- Moreover in general fs{j) £ ^j+s(j — 1) 
where j, s > 2. The basis monomials of the spaces Vy{^) in which we can express the func- 
tions fs{j) can be found, for example, in |25j . 

Proposition 1. If for each fixed j > 2 the equation ()2.14p with s = 2 and s' = 3, namely 

M2/3 (j) = A/3/2 (j) , (2.15) 

is satisfied, then there exist unique differential polynomials fs{j) Vs > 4 such that the flows 
MgU- = fs (j) commute for any s > 2 \24\ [6]. 

Hence among the relations (I2.14P only those with s = 2 and s' = 3 have to be tested. 

Proposition 2. The homogeneous equation MgU = has no solution u in the vector space V^, 
i.e. Ker(M^)nP^ = 0. 

Consequently the multiple scale expansion (|2.13p is secularity-free. This does not mean that, 
in solving equation ()2.13bp . we have to set to zero all the contributions to the solution coming 
from the homogeneous equation but only that part of it which is present in the forcing terms. 
Finally: 

Definition 4. If the relations (|2.14p are satisfied up to the index j, j > 2, we say that our 
equation is asymptotically C-integrable of degree j or Aj C-integrable. 

2.2.1 Integrability conditions for the Schrodinger hierarchy 

Here we present the conditions for the asymptotic C-integrability of order k or A/. C-integrability 
conditions with A; = 1, 2, 3. To simplify the notation, we will use for u- the concise form u{j). 

The Ai C-integrability condition is given by the absence of the coefficient p2 of the nonlinear 
term in the NLSE (fTT]) . 

The A2 integrability conditions are obtained choosing j = 2 in the compatibility condi- 
tions ()2.14p with s = 2 and s' = 3 as in ()2.15p . In this case we have that /2(2) S ^4(1) 
and /3(2) G "^5(1) where 7^4(1) contains 2 different differential monomials and 1^5(1) contains 
5 different differential monomials, so that /2(2) and /3(2) will be respectively identified by 2 
and 5 complex constants 

/2(2) = an„,(l)|tx(l)|2 + 6t2,„,(l)n(l)2, (2.16) 

/3(2) = a\uil)\\il) + /3|tx,„,(l)|\(l) +7n„,(l)2^(l) + 6u2nA'^)uilf + e|tx(l)|\,2ni(l). 

In this way, eliminating from equation ()2.15p the derivatives of u{l) with respect to the slow- 
times 7712 and m3 using the evolutions (I2.13ap with s = 2, 3 and equating term by term, we 
obtain that the A2 C-integrability conditions gives no constraints on the coefficients a and b 
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appearing in /2(2). The expression of the coefficients a, /3, 7, 6, e appearing in /a (2) in terms 
of a and b are 

a = 0, /? = -— , 7 = — --— , 6 = 0, £ = 7. 

i)2 ^-02 

The A3 C-integrability conditions are derived in a similar way setting j = 3 in equation (j2.15p . 
In this case we have that /2(3) G '^5(2) and /3(3) € 'Pq{2) where 7^5(2) contains 12 different 
differential monomials and 'Pq{2) contains 26 different differential monomials, so that /2(3) 
and /3(3) will be respectively identified by 12 and 26 complex constants 

/2(3)=Ti|n(l)|4«(l)+T2|«,„,(l)pu(l)+T3|n(l)|\,2ni(l)+r4n2ni(l)n(l)2 

+ T5Un, (l)'n(l) + Teun, (2) |t/(l) ^ + Tju^^ (2)n(l)2 + T8ui2)\{l) (2.17) 

+ T9|7i(2)|\(l) + Tiou(2)n „, (l)u(l) + TiM2)un^ (1)^(1) + ri2u(2)n ^^ (l)n(l), 
/3(3)=7i|u(l)|^n„,(l)+72|n(l)|\(l)2t2,„,(l)+73|t,(l)|2n3„,(l)+74tx(l)2n3„,(l) 

+ 75|'U,ni (1) P-W.ni (1) + l6U,2ni (l)«,ni (l)""!!) + l7U,2ni (l)^i,ni (l)""!!) 

+ 78^,2ni(l)n„i(l)72(l) +79^(1)1^^(2) +7io|n(l)|\(l)2^(2)+7nn„,(l)n(2)2 
+ 7i2W,ni (1) k(2) P + 713 k,ni (1) Pu(2) + 714 k(2) Pn(2) + 7i5U,„, (l)2n(2) 

+ 7l6|u(l)|2-u,2ni(2) +7l7u(l)^^,2ni(2) +718^(2)^,2^(1)^(1) 

+ 7i9w(2)n2ni (1)^(1) +720'u(2)n2ni(l)'u(l) +72iu(2)n„,(2)n(l) 
+ 722w(2)u „i (2)n(l) + 723^ m (2)n „i (l)n(l) + 724^711 (2)n „j (l)u(l) 

+ 725'",ni(2)n„i(l)'u(l) +726'U,ni(2)'u(2)n(l). 

Let us eliminate from equation ()2.15p with j = 3 the derivatives of u{l) with respect to the slow- 
times 7712 and m^ using the evolutions (I2.13ap with s = 2, 3 and the same derivatives of u(2) 
using the evolutions ()2.13bp with s = 2, 3. Equating the remaining terms term by term, the A3 
C-integrability conditions turn out to be: 

i _ - 1 

n = -7^- [b (rii - 2Tfi) + 0T7] , bTj = -{b- a) (m + no - Te) + arj, 

4ii2 2, 

arg = bTg = 0, org = brg = 0, aru = a (no - ni) + ftrg + arj, 

{b-a)n2 = {b-a)Tio. (2.18) 

Sometimes a and b turn out to be both real. In this case the conditions given in equations (|2.18p 
becomes: 

^1 = 7^ [^ (^11 - 2^6) + alj] , h = -7^ [b (Ru - 2Re) + aR^] , 

(6 - a) {Ru + i?io - ^6 - 2i?7) = 0, (6 - a) (In + ho -h- 2/7) = 0, 

(6 - a) i?8 = 0, (6 - a) /s = 0, (6 - a) iJg = 0, (6 - a) h = 0, 

a (i?i2 + -Rii - -Rio - Rj) = bRe, a (/12 + hi - ho - h) = bh, 

{b - a) {Ru - Rio) = 0, (6 - a) (/12 - Iio) = 0, (2.19) 

where ti = Ri-\- ili for i = 1, . . . , 12. The expressions of the 7j as functions of the n are: 

71 = jm (^■^e - 4i-B2n + &T12) , 72 = JH2 (^"^6 + CLTj) , 

3^^3X3 3l53T2 3iB3T4 

73 = 1^ — , 74 = U, 75 = ^^~' ')'6 = ^ — ' 



77 = 75, 78 = 


= 73 


SiB^Tg 

^''~ 2B, ' 




SiBsTQ 
^'' 2B2 ' 




SiBsTs 

^''= B, ' 




MB^Tj 
725 - „ , 
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5 — , 79 = 710 = 711 = 0, 

-D2 

713 = ^^ , 714 = U, 715 = — — , 

SiSgTio 

717 = 718 = U, 719 = TTB — ' 720 = 715, 

2i32 

722 = 712, 723 = 716 + 719, 724 = 7l3, 

726 = 0. 

The conditions given in equations (j2.18p . (j2.19p appear to be new. Their importance resides in 
the fact that a C-integrable equation must satisfy those conditions. 

3 Linearizability of the equations of the Burgers hierarchy 

Taking into account the results of the previous section we can carry out the multiple scale 
expansion of the equations of the Burgers hierarchy. To do so we substitute the definition (j2.1ip 
into equations (j2.ip , (j2.6p and write down the coefficients of the various harmonics 9 and of the 
various orders j of e. When we deal with the differential-difference equation (j2.ip . we have to 
make the substitutions am — ?> i, ami — )> tj. This transformation implies that in this case the 
corresponding coefficients p2 and B2 will turn out to be u-independent. In Appendix we present 
all relevant equations and here we just present their results. 

Proposition 3. The differential- difference equation (|2.ip of the Burgers hierarchy satisfies 
the Ai {and consequently also the A2) and also the A^ C -integrahility conditions. 



Proposition 4. The partial difference Burgers-like equation ()2.6p reduces for j = 3 and a = 1 

to a NLSE with a nonlinear complex coefficient p2 given by equation ()A.9cP . Thus the equation 
is neither S-integrable nor C-integrable. 

4 Conclusions 

In the present paper we have presented all the steps necessary to apply the perturbative multiple 
scale expansion to dispersive nonlinear differential-difference or partial difference equations 
which may be linearizable. These passages involve the representation of the lattice variables in 
terms of an infinite set of derivative with respect to the lattice index and the analysis of the 
higher order of the perturbation which give rise to a set of compatible higher order linear PDE's 
belonging to the hierarchy of the Schrodinger equation. The compatibility of these equations 
give rise to a linearizability test. We applied the so obtained test to the case of a differential- 
difference dispersive Burgers equation and its discretization. It turns out that the Burgers 
is linearizable (as it should be) but its discretization is neither S'-integrable nor C-integrable. 
So, effectively this procedure is able to distinguish between linearizable and non-linearizable 
equations. 

A Appendix 

Let us now start performing a multiple scale analysis of the partial difference equation ()2.6p . We 
present here the equations we get at the various orders of e and for the different harmonics 9. 

• Order e and 9 = 0: In this case the resulting equation is automatically satisfied. 
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• Order e and = 1: If one requires that tij^ / 0, one obtains the dispersion relation 
sin (k/i) a 



sin (w(t) 



h 



(A.l) 



• Order e^ and 6 = 0: We obtain the evolution 

0, 



Omi'^l jC/rii'^i 



(A.2) 



which implies that u\ depends on the variable p = hni + ami. 

Order e^ and 6=1: Taking into account the dispersion relation (jA.ip . we have 



9 u^^^ 

-'mi "1 



a 



cos (wo") 

(1) 



COS (k/i) 



(?„-, uv' 



-'rii ""l 



2sin^ 



2 /^^^^„(0).,{i) 



ul Ui 



0, 



which implies that u\ has the form 



u 



(1) 



5 (C, mj , j > 2) exp <j 5 / u^ ^ (p') dp 



,(0) 



(5 = 



2sin2(K;/i/2) 
[cos (k/i) — COS (wo-)] 



(A.3) 



, (A.4) 



where ^ = hni + ^^T^^crmi, (7 is an arbitrary function of its arguments going to zero as 



^ —7- ±00 and po) by a proper redefinition of g, can always be chosen to be a zero of u\ 



(0) 



,(0) 



as when p -^ ±00, u\ ^ so that there will exists at least one zero. 

Order e^ and 6 = 2: Taking into account the dispersion relation (jA.ip . we have 



u, 



(2) 



-inh 



l)h 



2 [cos {kK) — cos {uj(t)\ 



[u\ 



(1)12 



(A.5) 



Order e^ and 6 = 0: We have 






^mi^2 ' C/m Wo — C/m.^U 



(0) 

1712 "1 



2cr sin^ (k/i/2) [(9„j + 2/in[ 



(0)ll„(i)|2 



u 



(A.6) 



By equation (|A.2p , the term 9^2 i^i is a solution of the left hand side of equation ()A.6P , 
hence it is a secular term. As a consequence we have to require that 



dmiU2 — T9niU2 — — 2<7 siu^ (k/i/2) [8^ + 2hu\ 



(0)ll„(l)|2 



0-1712'^ I 



(0) 



h 
0. 



"1 I ' 



Solving equation (jAj) taking into account equation ()A.4p . we obtain 



u. 



(0) 



/ (/o, nij , j > 2) — M cos ((^(t) 






(A.7) 



where / is an arbitrary function of its arguments going to zero as p — ?• icxo and ^o is an 
arbitrary value of the variable ^. Let us restrict ourselves for simplicity to the case where 
there is no dependence at all on p. If one wants that the harmonic ul depends on £^ and 
not on p, from equation (IA.2P one has that d^ul = 0, so that u\ depends on the slow 
variables mj, j > 2 only. Similarly we have that d^f = 0. In this case, in order to satisfy 

the asymptotic conditions lim u\ = 0, 7 = 1, 2, one has to take u\ = / = (unless 

f— i-±oo 

we take the fully continuous limit /i — > 0, hni = xi, o" — )• 0, ami = h in which p ^ S,)- 
Equation (]A.7p then becomes 



u 



(0) 



-Mcos {uja) \ul 



{1)|2 



(A. 
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• Order e^ and 9 = 1: Taking into account the dispersion relation (jA.ip and the equa- 
tions uf^ = 0, (|X3]1 . (1X5]) . (|X8]1 . we have 

(1) (J cos (k/i) (1) (1) .„ ^2, (1) ,•„ L,(i)|2,,(i) ^A n«^ 

/lCOS(a;(T) ^ I 1 I 1 

/icr ih? — (T^) sin (k/i) 

^2 = ,r 2 - 2, M \21 V (A-9b) 

2 [(T^ sm (k/ij — /i^J cos \ojct) 

2ha sin^ (wcr/2) [2 cos (k/i/2) - cos (3k/i/2)] sin (k/i/2) 
^2 = 



[cos (k/i) — cos (wO")] COS ((JO") 

2i/icr sin^ {uia/2) sin (k/i/2) sin (k/i/2) 
[cos (k/i) — cos (wo")] cos (wcr) 



(A.9c) 



As a consequence of equation ()A.3P with n^ =0, the right hand side of equation ()A.9aP 
is secular. Hence we have to require that 

d) crcos (kK) „ d-i , . N 

^-1^2 - T y^dn^^ = 0' A. 10a 

n cos [uja) 

i9^, nj'^ = 5251^1'^ + P2 hS'^ I '^S'^ (A. 10b) 

Equation (jA.lOaP implies that Ug also depends on ^ while equation (jA.lObp is a nonintegrable 
nonlinear Schrodinger equation, as from the definition (IA.9cP we can see that p2 is a complex 
coefficient. So we can conclude that equation (j2.6p is not j4i-integrable. 

Let us perform the multiple scale reduction of the Burgers equation (12. ip . Equation (12. Ih can 
be always obtained as a semicontinuous limit of equation (j2.6p defining the slow times tj = anij, 
j > 1. In such a way we can use in the present calculation the results presented up above. 

• Order e and 6 = 0: In this case the resulting equation is automatically satisfied. 

• Order e and 9 = 1: Taking the semi continuous limit of equation (jA.ip . one obtains the 
dispersion relation 

^^_sin(K^_ (A.ll) 

h 

• Order e^ and 9 = 0: Taking the semi continuous limit of equation ()A.2p . one obtains 

dt,uf^-^dn,uf^ = 0, (A.12) 

which implies that u\ depends on the variable p = hni + ti. 

• Order e^ and a = 1: Taking the semi continuous limit of equation (jA.SP , one obtains 

a,^nW - S^iMa^^nW + 2sin2 {nh/2) nf ^n^ = 0, (A.13) 

which implies that u\ has the form 

u^i^ = 9^'^ (e, t„ J > 2) exp |- I" ^S°) (pO dp j , (A.14) 

where ^ = hui + cos {nh) t\, g^^' is an arbitrary function of its arguments going to zero as 
^ — >■ ±c» and po) by a proper redefinition of g, can always be chosen to be a zero ol u\ . 
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Order e^ and 9 = 2: Taking the semi continuous limit of equation (jA.Sp . one obtains 



u, 



(2) 



h 



-U 



(1)2 



Order e^ and 6 = 0: Taking the semi continuous limit of equation ()A.6p . one obtains 



(A.15) 



1 



,(0) 



,(0) 



5ti4 - -r5ni4 = -dt^ur - 2sin2 {Kh/2) [a„, + 2hu^^>] \uY^\. 



(0)lL,(l)|2 



h 



(A.16) 



By equation ()A.12p . the term 9*2 ^i i^ ^ solution of the left hand side of equation ()A.16p . 
hence it is a secular term. As a consequence we have to require that 



a, 



'ti^2 



(0) 



1 



5, 



ni^2 



(0) 



-2sin2(K/i/2) [5„i +2hu\ 



(0)ll.,(l)|2 



"1 I ' 



(A.17) 



dt.uf^ = 0. 



Solving equation (1A.17P taking into account equation (IA.14p . we obtain 



n, 



(0) 



,(i)|2 



(A.18) 



^2 - f {P^tj^J > 2) + h\u\ '\ , 

where / is an arbitrary function of its arguments going to zero p — ?> ±00. 

Order e^ and 9 = 1: For simplicity, from now on we require no dependence on 4j so that, 
in order to satisfy the asymptotic conditions, it necessarily follows that 



uf^ = f = 0. 
Taking the semi continuous limit of equations ()A.9aP , ()A.9bp , one obtains 

(1) cos(Kfe) (1) _ (1) 2, (1) R ^ h sin JKh) 



(A.19) 



(A.20) 



As a consequence of equation (|A.13p (with u\ =0), the right hand side of equation (IA.20p 
is secular. Hence we have to require that 



5.,n« - ^a„,n« = 0, 



(A.21a) 
(A.21b) 



Equation ()A.21aP implies that 1*2 depends also on ^ while, contrary to equation (jA.lObp . 
equation (IA.21bp now is a linear Schrodinger equation, reflecting the C-integrability of 
equation (|2.ip . 

Order e^ and a = 2: Taking into account the dispersion relation (lA.llh . the fact that 



u 



(0) 



and the equations (jA.lSp . (jA.lSp . we have 



uf^ = h 



1-e 



'.Kh 2 



(1) 



h 



4sm^{Kh/2) 



d^u\ 



(1) 



,(1) 



(A.22) 



^If 9puf ' / 0, dpf / 0, we have 

(1) • (2) / . • -^ i\ 

L ■J po 

with p arbitrary function of its arguments going to zero as ^ — > dzoo 






i f hih 



idtj^^ = B2dU^\ 



uf^+^Turdp'- rf{p')dp', 

J PO J PO 
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• Order e^ and 9=3: Taking into account the dispersion relation (|A.lip and equation ()A.15p . 
we obtain 

2 

,a)i3 



U 



(3) 



h 



1 -e* 



Kh 



H' 



(A.23) 



• Order e^ and 6 = 0: Taking into account equations (jA.lSp . (jA.lSp . (|A.19j) . (|A.21bp . we 
get 



5ti4°^ - j^dn.uf 



hd, 



h 



d, 



'«l^i 



(1)|2 



(1)„(-1)^„("1)„(1)^ 



2sin2(AcV2)(4^<^+4"^<0 



Solving equation ()A.24p . we obtain 



n. 



(0) 



r(p,tj,j >2) + /i(4^ui ' +u\ 'u\') 



(i).,(-i) , .,(-i).,(i)^ 



h^ 



-A\y>i 



(1)|2 



(A.24) 



(A.25) 



4sin [Kh/ 1) 

where r is an arbitrary function of its arguments going to zero as p — ?• ±00. As usually, 
if we don't want any dependence at all from p but only on ^, in order to satisfy the 

asymptotic conditions lim Ug = 0, we have to take 

g— >±oo 

r = (A.26) 

(unless we take the fully continuous limit). 

Order e^ and 9 = 1: Taking into account the dispersion relation (lA.lip and equations 
(|AT5l) . (lAlSJ) . (|AT9|) . (1X221) . (1X25]) . (1X261) . we get 

(1) COs(k/i) (1) (1) (1) 3 (1) 2, (1) 



5, 



^3 L '^"i ""3 

h"^ COS (k/i) 

6 



(A.27) 



As a consequence of the equations ()A.13p . ()A.21aP and as ul = 0, the right hand side of 
equation ()A.27p is secular, so that 



5.,nW-^a„,nW = 0, 



dt.u'i^ + iB2d'lu^2^ 



,(1) 



(1) 



B3dlu^^\ 



(A.28a) 
(A.28b) 



The first relation implies that Ug also depends on ^ while the second one, as a consequence 
of equation ()A.21bp . implies that the right hand side is secular, so that 



■fi (1) 



B2dl4^ 

l3.,(l) 



0, 



(A.29a) 
(A.29b) 



dt,u\> + Bzdlu\"=0. 

Equation ()A.29ap . as one can see from the definition (j2.16p . has a forcing term /2(2) with 
coefficients a = 6 = 0. 

Order e^ and 9 = 2: Taking into account the dispersion relation (jA.lip and the equa- 
tions ([Al3]) . (|AT5l) . (|AT8]) . ([Al9]) . dMH), ^EM, (1X23]) . we get 



n 



(2) 



oifth 



(1) (1)2 icos(KV2) 2 (1) 
7 til IUi I -r 7; — ; — 5-^ — ; — j:rr<^f ""i 



^"-1 I"-! 



/l^ 



8 sin^ (k/i/2) ' 
h 



,(1) 



4sin^(«;/i/2) 



a5("i"4") + rTi^(4"= + ^"4"). 



(A.30) 
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• Order e^ and 9 = 3: Taking into account the dispersion relation (|A.lip and the equa- 
tions (1X13]) . (|XT5]) . (IXT9]) . (r02]) . (1X231) . we get 

2 / ,., r,iiKh 



-(') - ' " u'^^^ 3n(^' + ^^^a,^(i) 



Order e^ and = 4: Taking into account the dispersion relation (jA.lip and the equa- 
tions (lATTSl) . dAiMI), we get 



ni^) = ( -^ f K^T- 



1-e 



IK/i 



Order e^ and 9 = 0: Taking into account equations fOS]) . (jXTSl) . (|AlT9]) . (jA.21bp . (|Al22]) . 
(1X261) . (TOSJ) . (1X281) . we get 

5^^ 40) _ ia„,nf = /ia^l - 2sin2 (k/i/2) (^I'^n^-^^ + nS-^)4^) + j^'^p) 

+ [4sin2 (Kh/2) - 1] y |4'Y + ^5c(^S'^4"'^ + ^^'^4'^) 
+ '-^cot{Kh/2)di:{u[~^'^d^u'i^ -u']^'^d^u["^'^)y (A.31) 

Solving equation ()A.3ip . we obtain 



u 



(°) = e {p,t„j > 2) + ^4^)4-^) + 4-^)4^) + i4^)p 



+ 



/l2 



a5(4'^4-'^ +4"')4'^) + [4sin2(KV2) - 1] /i|4'T} 



4:sm'^{Kh/2) 

ih^ cos {Kh/2) ( (~i)g (i) _ ^(l)g^^(-l)^ (j^ 32) 

8sin3(KV2) ^^ ' ^ ' ' ^ ' ^' ^ ^ ^ 

where is an arbitrary function of its arguments going to zero as p -^ ±cx). As usually, 
if we don't want any dependence from p but only on ^, in order to satisfy the asymptotic 
conditions lim u. = 0, we have to take 

6 = 

(unless we take the fully continuous limit). 

Order e^ and 9 = 1: Taking into account the dispersion relation (lA.llh and equations 
nj°^ = / = r = e = 0, together with the equations ([XTS]) . (lAlSJ) . (1X22]) . (|X23]) . (|X25]) . 
(1X30]) . (fX32l) . we get 

+ tB.d'A^^ + c [^S'^|5„,4')f + uA\dnA^f + 4'^'C4"'T > (A.33) 

. /i'^ sin (k/i) . /i [1 -|- cos (k/i) + 3i sin (k/i)] 

24 ^ gtKfc _ 1 

As a consequence of the equations ()A.13p . ()A.21ap . ()A.28aP and oiu\ = 0, the right hand 
side of equation (|A.33p is secular, so that 
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^4„,(l) 



W\r, „,(1)|2 , „,{-l)/o „,(1)\2 , „,(1)2;.2 „,(-!) 



The first relation tells us that u\ depends on ^ too while in the second one, as a con- 
sequence of equations ()A.21bp . ()A.29ap . the first our terms in the right hand side of 
equation ()A.34p are secular, so that 



,W^.-R./^2„(l) 



(1)1/^ „(1)|2^„(-1)^/^ ,,(1)^2 ^.,(1)2/^2 ,,(-1) 



9*3 4^^ + ^s^f-u^^^ = -duv^i^ + i549|'uS^^ (A.35b) 

As a consequence of equation ()A.29bp , the right hand side of equation ()A.35bp is secular 
so that 

dp,u^2^ + Bsdlu^^^ = 0, dt^u^^^ - iBid^v^l^ = 0. 

Equation ()A.35ap . as one can see from the definition (|2.17p and taking into account 
that a = 6 = 0, has a forcing term /2(3) that respects all the A^ C-integrability con- 
ditions (f2J8]l . 
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